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Supernova Neutrino Oscillations * 

H. Nunokawa a 

a Instituto de Ffsica Gleb Wataghin, Universidade Estadual de Campinas, UNICAMP 
13083-970 Campinas SP, Brazil 

We discuss some possible influence of neutrino oscillations on physics of supernova mainly focusing on the 
observations of u £ , and present some analysis of SN1987A data in the light of three neutrino mixing scheme. 



1. Introduction 

Sucessful observation of neutrino burst from 
the supernova SN1987A, which is believed to be 
the collapse driven one, in the Large Magellanic 
Cloud by Kamiokande-II Q and 1MB g detec- 
tors have confirmed the basic picture of supernova 
explosion (see Refs. |||4| for a review). 

Such collapse driven supernova can be a good 
laboratory to test various unknown neutrino 
properties not only mass and flavor mixing but 
also other properties such as magnetic moment, 
life time and some other non-standard interac- 
tions j|. 

It has been known that neutrino flavor con- 
version by the MSW effect @ inside supernova 
could cause some significant influence on super- 
nova physics |§J^-§). In this talk we review the 
possible influence of neutrino flavor conversion on 
supernova physics and present some analysis of 
the SN1987A data in light of three flavor mixing 
scheme of neutrinos, based on our recent work Q. 

2. Standard picture of supernova neutrino 
emission 

Here we summarize the basic understanding of 
neutrinos from supernova (SN) Band their prop- 
erties inside neutrinosphere JlOj — 11 2|] - 

A type-II supernova occurs when a massive star 
(M > 8Mq) has reached the last stage of its life. 
Almost all (~ 99%) of the gravitational binding 
energy of the final neutron star (about ~ 10 53 erg) 
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is radiated away in form of all flavors of neutri- 
nos. The individual total neutrino luminosities in 
supernovae are approximately the same for all fla- 
vors but the individual neutrino energy distribu- 

This is because dif- 



tions are very different! 10 



li- 



ferent flavor of neutrinos interact differently with 
the star material, as implied by the following re- 
actions, 



v e + n 

i> e +P 
v + N 



,+ 



V 

n + e 1 
v + N, (N=p,n). 



(1) 
(2) 
(3) 



Since the cross sections of the charged-current 
reaction is larger than that of the neutral-current 
one and there are more neutrons than protons, 
the i/ e 's have the largest interaction rates with 
the matter and hence thermally decouple at the 
lowest temperature. On the other hand, v T (p) and 
v T (u) ' s l ac k the the charged-current absorption re- 
actions on the free nucleons inside the neutron 
star and hence thermally decouple at the highest 
temperature. 

As a result, the average neutrino energies sat- 
isfy the following hierarchy: 



(E Ve )<(E B .)<{E VTt Jn{E p . ) ). 



(4) 



Typically, the average supernova neutrino ener- 
gies are, (E„ e ) w 11 - 12 MeV, {E De ) a 14- 
17MeV, (E VrM ) w (Eb r00 ) « 24-27 MeV. Such 
energy hierarchy is of crucial importance for our 
discussion as we will see later. 

The shape of the energy spectra of various fla- 
vors of neutrinos can be described by a " pinched" 
Fermi-Dirac distribution |Q. The pinched form 
can be parametrized by introducing an effective 
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"chemical potential", r\ as follows, 
F{E V ) oc 



(5) 



1 + exp[-E v /T - rj\ ' 

where T is the temperature. 

There is no physical, significant distinction be- 
tween and v T and their antiparticles in neu- 
trinosphere. It is because and are not en- 
ergetic enough to produce muons by the charged 
current interactions, and the neutral current cross 
sections of v and v are similar in magnitude. 
Therefore, following Ref. |9) we collectively de- 
note them as "heavy falvor neutrinos", Vheavyi hi 
this talk even if they are not heavy (as in the case 
of nverted mass hierarchy we will discuss later). 

3. Possible influence of neutrino oscillation 
on supernova physics 

Here we briefly review possible effects of neu- 
trino oscillation on supernova physics, which were 
discussed in number of previous works. 

First let us stress some characteristic features 
of neutrino oscillations in supernova. Because 
of extremely high matter density around neutri- 
nosphere, MSW resonant neutrino conversion ]j| 
can occur for the mass squared difference much 
larger than the ones relevant for solar or atmo- 
spheric neutrinos. Because of this, in general, one 
must deal with multiple resonances in supernova 
(in contrast to the case of solar neutrinos) in the 
case of three (or more) flavor mixing scheme. 

Another important point is that because super- 
nova emit all flavors of neutrinos, if v a — > vp or 
v a — > vp (a, = e, heavy) occurs, vp — > v a or 
Dp — ► D a also must occur simultaneously, and the 
net effect in the case of large oscillation probabil- 
ity is the exchange of the energy spectrum of v a 
and vp or D a and Pp. 

3.1. Shock reheating 

Neutrino conversion could have significant im- 
pact on shock re-heating in the delayed explo- 
sion scenario JL5| ]. If the conversion between v e 
and Vheavy occurs in the region between the neu- 
trinosphere and the stalled shock this can help 
the explosion 0. Due to the conversion the en- 
ergy spectra of v e and Vheavy can be swapped and 
hence v e would have larger average energy leading 



to a larger energy deposition by reactions in eqs. 
(|l|) and (||) so that the stalled shock would be re- 
energized. In order for this effect to be operative, 
Am 2 > 100 eV 2 is required @. 

3.2. Heavy elements nucleosynthesis 

In nature, supernova is considered to be one 
of the most promising site to create neutron rich 
heavy elements [|l6| . To have successful r-process 
the site must be neutron rich, i.e. Y e < 0.5 where 
Y e is number of electron per baryon. The Y e value 
is mainly determined by the competition between 
the two absorption reactions in eqs. (|IJ) and (0). 
In the standard supernova model the latter pro- 
cess is favoured due to the higher average energy 
of D e which guarantees the neutron richness. If 
the neutrino oscillations in v e — Vheavy channel do 
occur between the neutrinosphere and the region 
relevant for r-process the site can be driven to 
proton-rich due to the reaction (||) and therefore, 
r-process could be prevented ||7 Using this ar- 
gument, mixing parameters in the region Am 2 > 
few eV 2 and sin 2 29 > 10~ 5 - 10~ 4 can be ex- 
cluded [§. 

3.3. Influence for v e induced events 

If v e — Vheavy oscillation occurs, some enhance- 
ment of forward peaking elastic scattering events 
at high energies which should be observable in 
water Cherenkov detectors flffl, as well as en- 
hanced oxygen-induced events due to a steep rise 
of the cross section at energies higher than > 
30 McV pgj] , which could be separated from the 
dominant isotropic i7 e absorption events due to a 
moderate backward peaking of the events |19|. 



3.4. v e signal 

Last argument is the effect of neutrino oscilla- 
tion on v e signal in the terrestrial detector, which 
will be discussed more in detail in the follow- 
ing sections. If large oscillation between v e and 
Pheavy occurs, the v e spectrum gets harder and 
can affect significantly the observation of super- 
nova neutrinos, because the cross section for the 
v e p absorbed charged current reaction is much 
larger than that of the v e e~ elastic one. Apply- 
ing this argument, vacuum oscillation and large 
mixing angle MSW solutions to the solar neutrino 
problem are disfavored by the observed SN1987A 
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neutrino data at Kamiokande and 1MB ||] de- 
tectors ||^-|2^]. We will discuss in sec. 5 that 
how this conclusion could be modified or to be 
interpreted in the context of three neutrino mix- 
ing scheme. 

4. Neutrino conversion in three flavor 
scheme 

We now discuss the neutrino flavor conversion 
in supernova (SN), in the three flavor mixing 
scheme. To be most conservative, here we assume 
neutrino mixing only among three active neutrino 
flavors, i.e., among v ei v T and among their 
anti-particles, and do not consider the conversion 
into sterile neutrinos. 

Here, we assume oscillation interpretation of 
the atmospheric |2^] and solar neutrin o |24| l data 
and do not consider the LSND result p5|. Some 
possible implications of the mixing scheme which 
can explain the LSND result have been discussed 
in Refs.||-|§. 

We note that two Am 2 which can explain at- 
mospheric as well as solar neutrinos are irrelevant 
for the shock reheating and heavy elements nucle- 
osynthesis (see sees. 3.1 and 3.2) and therefor, we 
will discuss the influence only on v e and v e sig- 
nals, mainly focusing on the latter. 

The first very important question is if the neu- 
trino mass spectrum adopts the normal or in- 
verted mass hierarchies, because it determines 
wheather there is a resonant conversion in the 
antineutrino channel or not. These two mass hi- 
erarchies are illustrated schematically in Fig. 1, 
in which we assume that smaller mass squared 
difference is in the range relevant to explain the 
solar neutrino data and the larger one for atmo- 
spheric neutrino data. 

In these mass hierarchies, the three neutrino 
and three antineutrino eigenstates have two level 
crossings, first at higher (H) density and the sec- 
ond at lower (L) density, inside SN as schemati- 
cally illustrated in Fig. 2. If the mass hierarchies 
is of normal (inverted) type, the H level crossing 
is in the neutrino (antineutrino) channel. 

The second important question is that if the 
neutrino (or antineutrino) flavor conversion in SN 
at H level crossing is adiabatic or not. If it is very 



(a) Normal 



(b) Inverted 



Am. 
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Fig. 1: Mass hierarchy schemes we consider. 



adiabatic, then the physical properties of neutrino 
conversion is simply v e - v heavy or v e - v h eav V ex- 
change in the case of normal or inverted mass 
hierarchy. It should be emphasized that this fea- 
ture holds irrespective of the possible complexity 
of the solar neutrino conversion which governs 
the L resonance. These key features have been 
pointed out in our earlier paper, Ref. |]l7|| in the 
case of normal mass hierarchy. 

If the mass hierarchy is inverted and H reso- 
nance is adiabatic, then v e spectrum gets harder. 
Since the P e -induced charged current reaction is 
dominant in water Cherenkov detector, one can 
severely constrain the scenario of inverted mass 
hierarchy by utilizing this feature of neutrino fla- 
vor transformation in SN. When the next super- 
nova event comes it can be used to make clear 
judgement on whether the inverted mass hierar- 
chy is realized in nature, a completely indepen- 
dent information from those that will be obtained 
by the future long-baseline neutrino oscillation 
experiments |29-3l|, or by the neutrinoless dou- 
ble beta decay experiments ]33| for the case of 
Majorana neutrinos. 

The adiabaticity of the H resonance is guaran- 
teed if the following adiabaticity parameter 7 is 
significantly larger than unity at the resonance 
point: 



Am 2 sin 2 29 
2E cos 26» 



d In TV 



dr 



Am 2 \ 1 ~ 1/n sin 2 29 



2E 



(cos26>) 1+1/ ™ n 



V2GFPoYe 



1/ri 

(6) 
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(a) Normal 



(b) Inverted 




Fig. 2: The schematic level crossing diagram for the case of (a) normal and (b) inverted mass hierarchies 
considered in this work. The circles with the symbol H and L correspond to resonance which occur at higher 
and lower density, respectively. Adopted from Ref. M 



Here, we assumed that the density profile of the 
relevant region of the star can be described as 
p(r) = po(r /r )" n to obtain the second line in 
the above equation, where r = 6.96 x 10 10 cm 
denotes the solar radius. With the choice n = 3 
and p — 0.1 g/cc p3|, we get, 



7 ~ 0.63x 



'13 



10" 



A; 



10" 3 eV^ 



2/3 



E 



20 MeV 



-2/3 



,(7) 



for the small value of #13. Since the conversion 
probability P is approximately given by Pr — 
exp[— §7], sin 2 0i3 > a few x 10~ 4 assures adia- 
baticity in a good accuracy. 

For a recent complehensive treatment of neu- 
trino flavor conversion in SN in the framework of 
three-flavor mixing, see Ref. |p3. 

5. Analaysis of the supernova SN1987A 
neutrino data in the light of three neu- 
trino scheme 

While waiting for the next galactic SN, let us 
perform an analysis of the data of neutrinos from 
SN1987A to gain a hint to the problem of the 
mass and mixing pattern which we want to solve. 

We repeat the similar analyses performed in 
Refs. but in the context of three-flavor 

mixing scheme of neutrinos, which is essential for 
the SN neutrinos. In due course, we will try to 



clarify how conclusions obtained in earlier works 
are to be interpreted, or to be conditioned in the 
three-flavor framework. 

We follow the statistical analysis of the 
SN1987A data performed by Jegerlehner, Neu- 
big and Raffelt [glj who employed the method 
of maximum likelihood. Using the same Likeli- 
hood function which can be found in Ref. , we 
perform a fit of the observed data to the two pa- 
rameters, the binding energy of the neutron star, 
-Efc, and the u e temperature, T Pc . 

First we show the result of our analysis with- 
out assuming neutrino oscillation. In Fig. 3 
we show the contours of constant likelihood for 
Kamiokande and 1MB data, as well as the com- 
bined one. 

For simplicity, as in Ref. |21j) , we set the "ef- 
fective" chemical potential equal to zero in the 
neutrino distribution functions because we be- 
lieve that our results would not depend much even 
if we introduce some non-zero chemical potential. 
We note that our result is in very good agreement 
with the one obtained in Ref. 

We note that in Fig. 3, the agreement between 
Kamionde and 1MB data is not so significant 
though these data are consistent. The possibil- 
ity to reconcile the difference between Kamionde 
and 1MB data by utilizing the earth matter effect 
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has been recently considered in Ref. |35(]. Here, 
we do not enter into this point and simply try to 
do the combined fit as in Ref. Ell. 
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Fig. 3: Contours of constant likelihood corresponding 
to 68.3, 90 and 95.4 % C.L. region in E b - Tfc c for 
Kamiokande and 1MB data. Best fitted points are 
indicated by open circles. Think dashed curve shows 
the combined 95.4 % C.L. contour and the solid circle 
indicates the combined best fit. 

5.1. Case I 

Let us first consider the case where the observ- 
able effect is most significant: the mass hierar- 
chy is inverted and H resonance is very adiabatic, 
which can be realized if sin 2 6*13 > 10 -4 . fn this 
case there is a efficient conversion in v channel at 
H resonance. 

We draw in Fig. 4 equal likelihood contours 
as a function of the heavy to light tempera- 
ture ratio t on the space spanned by v e tem- 
perature and total neutrino luminosity by giving 
the neutrino events from SN1987A observed by 
Kamiokande [Q and fMB Q detectors. We char- 
acterize the difference in energy by the tempera- 
ture ratios of v e and v e to Vheavv as 



T, 



T = 



Vheavy 



~r (8) 

According to the simulation of supernova dynam- 
ics which is carried out in Ref. JTo|-[l^|, typically 
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t ~ 1.4-2.0. 
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Fig. 4: Contours of constant likelihood which corre- 
spond to 95.4 % confidence regions for the inverted 
mass hierarchy under the assumption of adiabatic H 
resonance. From left to right, r = T Bhe /T Be = 
T Vheavy /T Bt = 2,1.8,1.6,1.4,1.2 and 1.0. Best-fit 
points for T Ve and Eb are also shown by the open 
circles. The parameter rj parametrizes the departure 
from the equipartition of energy, L Vheav = Lp heav = 
i]L Vf , — rjL Vc , and the dotted lines (with best fit in- 
dicated by open squares) and the dashed lines (with 
best fit indicated by stars) are for the cases r\ — 0.7 
and 1.3, respectively. Theoretical predictions from su- 
pernova models are indicated by the shadowed box. 
Adopted from Ref. || 



In addition to it we introduce an extra parame- 
ter 77 defined by L Vheavy = L„ heavy = j]L Vp = r)L^ 
which describe the departure from equipartition 
of energies to three neutrino species and exam- 
ine the sensitivity of our conclusion against the 
change in the SN neutrino spectrum. 

At r = 1, that is at equal v e and v e temper- 
atures, the 95 % likelihood contour marginally 
overlaps with the theoretical expectation jl2j rep- 
resented by the shadowed box in Fig. 4. When 
the temperature ratio r is varied from unity to 2 
the likelihood contour moves to the left, indicat- 
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ing less and less consistency, as r increases, be- 
tween the standard theoretical expectation and 
the observed feature of the neutrino events after 
the MSW effect in SN is taken into account. This 
is simply because the observed energy spectrum 
of v e must be interpreted as that of the original 
one of Dheavy , in the presence of the MSW effect in 
the anti-neutrino channel, which implies that the 
original v e temperature must be lower by a fac- 
tor t than the observed one, leading to stronger 
inconsistency at larger r. 

The solid lines in Fig. 4 are for the case of 
equipartition of energy into three flavors, 77 = 1, 
whereas the dotted and the dashed lines are for 
r\ = 0.7 and 1.3, respectively. We observe that 
our result is very insensitive against the change 
in r\. 

We conclude that if the temperature ratio r is 
in the range 1.4-2.0 as the SN simulations indi- 
cate, the inverted hierarchy of neutrino masses 
is disfavored by the neutrino data of SN1987A 
unless the H resonance is nonadiabatic. We 
note that similar conclusion is obtained also in 
Ref. f|. 

5.2. Case II 

Let us next consider the case where mass hi- 
erarchy is normal or inverted but H resonance is 
very non-adiabatic, which correspond to the case 
with very small #13, sin 2 #13 <C 10~ 4 . In this case, 
there is no (or no significant) adiabatic conversion 
in v channel, and therefore, as far as v e signal is 
concerned, the problem is essentially reduced to 
that in the two flavor scheme which was studied 
in Rcfs. |0|j2|]. 

Conclusions can be summarized as follows: If 
the small mixing angle (SMA) MSW is the solu- 
tion to the solar neutrino problem there is only 
a minor effect because neither the vacuum oscil- 
lation nor the earth matter effects are effective 
because of small &n. If the large mixing angle 
(LMA) or low Am 2 (LOW) MSW or vacuum os- 
cillation (VO) is the solution, we have a potential 
trouble because a good fraction of v e is trans- 
formed into Pheavy an d vice versa. 

We have repeated the same analysis as done 
in Ref. pH and obtained very similar results. We 
conclude that VO and LOW solutions are more 



disfavored than the LMA MSW solution mainly 
because of the absence of earth matter effect for 
the former solutions as discussed in Refs. [po| , pT| . 
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Fig. 5: Contours of constant likelihood corresponding 
to 95.4 % C.L. for the large mixing angle MSW solu- 
tion (a) without and (b) with earth matter effect. We 
have taken mixing parameters as Am 2 = 3 x 10~ 5 eV 2 
and sin 2 26 = 0.8 for LMA MSW solution. Adopted 
from Ref. |] 

Let us try to demonstrate explicitly this point. 
We employ a particular set of parameters of the 
LMA MSW solution and compare the behavior 
of the likelihood contours with and without earth 
matter effect. We show our results in Fig. 5. For 
simplicity, we set 6> 13 = but our result does not 
change much as long as the parameter is under the 
CHOOZ bound |3(|. We see that inconsistency 
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between the fitted data and theoretical prediction 
is somewhat weakened when earth matter effect 
is included. 

5.3. Case III 

Finally, let us consider the case where mass hi- 
erarchy is inverted but H resonance is moderately 
adiabatic, which means that Pjj is not so close to 
zero or to unity, which can correspond to the case 
with sin 2 13 - 1CT 4 - 1CT 5 . 

In this case, v e -Vh eavy transformation occurs 
with the probability 1 — Pjj, and it would imply 
the similar but milder effect than that we have 
obtained with a good adiabaticity of the H res- 
onance, shown in Fig. 4. If the next galactic 
supernova is detected by Superkamiokande, then 
we will be able to discriminate the moderately 
nonadiabatic case from the adiabatic one. 

6. Summary 

We have discussed the possible influence of neu- 
trino oscillation for supernova physics, in partic- 
ular for observation of D e spectrum. We stress 
that the mass spectrum can be tested by v e sig- 
nal from supcrnovae, if #13 is not very small. We 
performed some analysis of SN1987A data in the 
context of three flavor mixing and conclude that 
inverted mass hierarchy is disfavored by the data 
unless #13 is very small, sin 2 6*13 < 10~ 4 . We hope 
that future galactic supernova may provide more 
clear information. 
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